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THERMOCHEMICAL CYCLE ANALYSIS USING LINKED 


CECS72 AND HYDRGN COMPUTER PROGRAMS 
by Leo F. Donovan 
Lewis Research Center 

SUMMARY 

The Lewis Research Caster chemical equilibrium computer program was 
coupled to the University of Kentucky thermochemical cycle analysis computer 
program. Input to the combined program is the same as input to the thermo- 
chemical cycle analysis program except that the extents of the reactions need 
not be specified. 

The combined program is designed to be run interactively from a computer 
time- sharing terminal. This mode of operation allows correction or modification 
of the cycle to take place during cycle aralysis. 

A group of 13 thermochemical cycles was used to test the combined program. 
Six cycles were found defective because one of the steps in each cycle was not 
feasible. Results of the analyses of the remaining seven cycles were in general 
less favorable than reported values. 


INTRODUCTION 

Closed thermochemical cycles have been proposed as a means of producing 
hydrogen and oxygen from water using a heat source, such as a nuclear reactor, 
for the thermal energy. Thermochemical cycles are currently of interest be- 
cause they may provide i more efficient alternative to electrolysis. 

Many new thermochemical cycles have been proposed in the recent literature 
(e. g. , refs. 1 to 3). A number of ways for deciding if these cycles warrant fur- 
ther study have been developed. They do not require the detailed calculations and 
experimental data needed for an engineering design but are helpful as a prelimi- 
nary screening to identify possibly interesting cycles and to eliminate those cycles 
with less potential. Some methods use only heats and free energies of formation. 
More involved methods require knowledge of other thermodynamic data as well as 
extents of reaction. Unfortunately, experimental values for extents of all reac- 
tions in a cycle are seldom known and estimated values must be used. Most cy- 
cles that have been analyzed for efficiency in this manner have supposedly used 
equilibrium extents of reaction. 
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The computer program HYDRGN was developed by Funk, et al. (ref. 4) for 
thermochemical cycle analysis. The program takes as input for each step in the 
cycle the chemical reaction occurring, the extent of the reaction, mole ratios of 
reactants, and temperature and pressure. The cycle evaluation program uses 
this information and appropriate thermodynamic data to calculate the work of 
separation of gaseous products, net heat required after as much heat exchange 
as possible between streams in the cycle, and finally, two figures of merit or 
efficiencies for the entire cycle. Results of calculations using this program have 
been reported in the literature (refs. 2, 5, and 6). 

A sequence of computer programs for calculating thermodynamic equilibria 
has evolved over the years at the NASA Lewis Research Center. The latest ver- 
sion of this program (similar to the program documented in ref. 7), referred to 
as CECS72, has the ability to handle solid and liquid products but not aqueous 
solutions. Thermodynamic data included with the progra. l are mostly taken from 
the JANAF tables (ref. 8). 

Current interest in energy systems led us to the decision to link the existing 
chemical equilibrium computer program to the University of Kentucky 1 s thermo- 
chemical cycle analysis computer program. This linking will allow checking of 
published cycles and facilitate evaluation of new cycles. For those steps in a 
cycle for which experimental values for extents of reaction are available, the 
equilibrium calculations can be bypassed. 

Since both HYDRGN and CECS72 are occasionally updated, it was decided 
to link them in such a way that the integrity of neither program would be affected. 
Any information passed between them would be via data sets that are written by 
one program and read by the other. In this way updates of either program could 
be incorporated with a minimum of trouble, but the original linking would be 
more difficult. 


LINKING PROGRAM 

The original plan was for a rather straightforward linking of CECS72 and 
HYDRGN. It was only during evaluation of the first test cases that it became 
apparent that the user must easily be able to modify the cycle. Three kinds of 
problems dictated the final design. The first of these was computational; the 
other two were in cycle design. 

Thermochemical cycles involve at least one solid or liquid in almost every 
reaction. In order for CECS72 to calculate equilibria with condensed species 
present, it is sometimes necessary to specify which of several possible con- 
densed species the program should consider first (using the " insert" feature). 
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We have found it is better to consider all condensed products of a reaction first. 
However, to accommodate the cases in which this is not successful, the user is 
allowed to specify another choice of insert (including no insert at all). 

At times products were encountered that were not anticipated - or were ig- 
nored - by a cycle designer. In some cases these extra products could be con- 
sumed as reactants in other steps. In other cases, the reactions were very 
temperature dependent; a 100 K change could alter the equilibrium composition 
significantly. It was found that the temperature of some steps needed to be 
changed in order to obtain the desired products. 

Figure l shows a sketch of the flow of the program. The symbols used in 
the linking program are contained in appendix A and a listing of the program is 
given in appendix B. The program was designed to be run interactively from a 
time- sharing terminal. The user is notified whenever a difficulty arises and a 
course of action to be taken is suggested. 

The linked programs, which are written in FORTRAN IV, run on an IBM 
360/67 computer operating under TSS. Extensive use is made of system facili- 
ties, especially the command procedure facility (ref. 9) and the research editor 
(ref. 10). All changes to data sets are made with the research editor. 

The linking program (which is actually two main programs), CECS72, and 
HYDRGN are executed from a procdef (CONTROL) requiring two parameters. 
The first parameter is the name of the input data set (in the form required by 
HYDRGN except that the extent of reaction need not be given). If the extent of 
reaction is known experimentally CECS72 can be bypassed for that reaction by 
specifying the second parameter. If this parameter, which acts as a two-way 
switch, is not blank program flow branches to a point immediately after the 
statement which executes CECS72. If the parameter is blank control passes to 
the next statement in the procdef. . 

The first linking program (CONI) read; CONREA) the inptit data set, re- 
formats (CONREF) each of the species names to be consistent with CECS72, and 
prepares (CONPRE) the input data set for CECS72. If any equation is not bal- 
anced, a message is typed at the terminal stating the problem and advising the 
user to correct the input data set and rerun the problem. Corrected data sets 
are checked until no more errors are detected, 

CECS72 is run to determine the equilibrium composition for those reactions 
in which the extent of reaction was not specified. The second linking program 
(CON2) is next run. After again reading and reformatting the input data set, the 
CECS72 output data is read (CONS72) and checked for products not originally 
specified (CONEXT), the extents of reaction are calculated (CONE PS), and the 
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original input data set is modified to incorporate the equilibrium extents of reac- 
tion (CONMOD). 

A final check for balanced equations is made. Also, the user is notified 
(default values are given here in parenthesis) if the concentration of unexpected 
products is too large (>6 mole percent) or if the extent of reaction is too small 
'^0. 01). ff so, the keyboard is unlocked and the user is invited to change the 
input data set and rerun the problem. 

ff the previous checkpoint is successfully traversed CECS72 output is 
checked for nonconvergence, ff the equilibrium calculation for a reaction was 
not successful the user is advised to change inserts and rerun with the second 
parameter not defaulted. 

Finally, HYDRGN is run to perform the cycle analysis and print the results. 
Also, printed in a separate listing for reference, are the warning messages that 
were printed at the terminal, a summary of CECS72 output, and a table of ex- 
tents of reaction. 


CYCLE EVALUATION 

In a recent paper (ref. 2) Funk combined 26 reactions in different ways to 
create 13 new thermochemical cycles which were evaluated using HYDRGN. The 
high temperatures required by these cycles are beyond current nuclear reactor 
capability (process heat temperature of about 1100 K), but may be available in 
the future. These cycles were reevaluated to see how well the linked programs 
worked. 

In the course of reevaluating these cycles a number of changes needed to be 
made to the reactions given in reference 2. 

1. hi some high temperature steps sulfur is given as the gaseous monatomic 
product S although equilibrium calculations indicate that Sg is the correct form. 
Substitution of S 2 for S requires the introduction of another step to condense S 2 
to S(l) in order to balance the cycle. This step can be carried out at 718 K, the 
boiling point of sulfur. 

2. Some cycles contain a low temperature step in which S(f ) appears as a 
product and a high temperature step involving S as a reactant. If one attempts 
to balance the cycle by introducing a step to vaporize S(l) at 850 K, the temper- 
ature recommended in reference 5, equilibrium calculations indicate that approxi- 
mately equimolar amounts of S 2 and S g are formed. Fortunately, however, the 
S 2 and Sg can be used as reactants in place of S in the high temperature step. 

Z. In some cases, the reaction temperature needed to be changed in order to 
reduce the amount of side products to an acceptable level. 
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4. One reaction involved CdS0 4 (l) decomposition at 1300 K. Since we have 
no data on liquid CdS0 4 , we substituted CdS0 4 (s) at its melting point, viz. , 

1273 K. 

The cycles as originally reported and as revised are shown in table I. 

Two reactions were unacceptable, thereby eliminating six cycles. The 
equilibrium pressure for 

1000 K 

Bi 2 S 3 (s) «... _ : 2 Bi(s) + 1. 5 S 2 

was only 0. 0028 atmosphere and 10 percent BiS was formed. For 

1273 K 

CdS0 4 (s) CdO(s) + S0 2 + 0. 5 0 2 

The extent of reaction was negligible («0. 01). 

A comparison of the results with those reported (ref. 2) for the seven cycles 
that were successfully calculated is given in table II. Molar recycle ratio for a 
step is the number of moles of gas that must be recycled at that step per mole of 
hydrogen produced by the cycle. The revised maximum molar recycle ratio is 
generally larger than that reported. Work of separation is calculated for gases 
only and assumes 100 percent separation efficiency. The reported and revised 
values are generally in agreement. Heat required after matching is the external 
heat input to the cycle. The reported and revised values agree quite well. 

The rj f s are estimates of thermal efficiency, assuming an efficiency of 
30 percent in converting heat to work. The first estimate, rjp assumes total 
recovery of all exothermic heat and 100 percent efficiency for the work of sepa- 
ration. The reported and revised rjj's both average 0. 53. In calculating the 
second estimate, r) 2 , it is assumed that the work of separation is 50 percent effi- 
cient and that work is generated from the exothermic heat (above 500 K) remain- 
ing after heat exchange matching. The revised n 2 ’s are lower than the reported 
ones, averaging 0. 31 as opposed to 0. 36. 

Some further information can be extracted by examining the detailed calcula- 
tions. Table III shows the molar recycle ratio, work of separation, and extent 
of reaction, c, for each step in the cycles successfully calculated. Some cycles 
can be improved significantly if the extent of reaction can be increased for cer- 
tain steps. For example, the reaction 

1300 K 

h 2 s^::z^h 2 + o.5S 2 
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which occurs as step 3 of cycle 1 and step 2 in cycles 5, 7, and 10 consumes 
50 to 60 percent of the 17 to 24 keel of separation work. In addition, this reac- 
tion accounts for one- third to one- half of the total recycle ratio. The reaction 

500 K 

0. 25 BaS(s) + H z O 7 » 0. 25 BaS0 4 (s) + H 2 

which occurs as step 1 of cycles 3 and 8, is responsible for almost die entire 
recycle ratio. Also, this reaction uses one-third to one-half of the work of 
separation of 7 to 11 keal. 

In order to determine the influence of the extents of reaction on cycle effi- 
ciency, the cycles that were successfully calculated using equilibrium values 
for extents of reaction were recalculated assuming complete conversion of re- 
actants to products. The results are given in table IV. Examination of the 
table demonstrates that erroneous conclusions about cycle efficiency can be 
drawn unless reliable values for extents of reaction are known. 

A simplified estimate of cycle efficiency that has been proposed (personal 
communication from M. G. Bowman) can be written as 



* H H 2 0, 298K 


AH 298K + 


pos 

Y AGp/0. 
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where the first summation is over all endothermic reactions at 298 K and the 
second summation is over all reactions with positive free energy changes at the 
temperature of the reaction. A comparison of this estimate with values calcu- 
lated with the HYDRGN program using equilibrium values for extents of reaction 
is given in table V. This comparison indicates that efficiencies obtained from 
this estimate are probably overly optimistic. 


CONCLUDING REMARKS 

The Lewis Research Center's chemical equilibrium computer program 
CECS72 was coupled to the University of Kentucky s thermochemical cycle 
analysis computer program HYDRGN. Input to the combined program is the 
same as input to HYDRGN except that the extents of reaction need not be speci- 
fied. After CECS72 calculates equilibrium compositions for each reaction, the 
extents of reaction are calculated by the linking program and inserted into the 
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input data set for HYDRGN. If experimental values are known for some or all 
of the reactions, they can be specified in the input and the equilibrium calcula- 
tions will be bypassed. 

The combined program is designed to be run interactively from a computer 
time- sharing terminal. If the program encounters errors in the input or cycle 
design flaws, the user is notified at the terminal and invited to correct or modify 
the offending data set. Reactants and products can be added to or deleted from 
existing reactions, new reactions can be added, or temperatures can be changed 
to modify a defective cycle into a potentially useful one. 

A group of 13 thermochemicai cycles constructed from 26 chemical reac- 
tions were used to test the combined program. Equilibrium calculations showed 
that one reaction had an unacceptably low equilibrium partial pressure and one 
reaction had a negligible extent of reaction; six cycles were thus eliminated. 
Results for the remaining seven cycles were generally less favorable than re- 
ported values. The steps that would have the greatest impact in improving the 
cycles were pointed out. 
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APPENDIX A 

SYMBOLS USED IN COMPUTER PROGRAM 

F4CLC, F4TRT, and F4MVC are FORTRAN implementations of the assem- 
bly language commands CLC (compare logical character), TRT (translate and 
test), and MVC (move character). 

The names of variables are the same as those used in HYERGN, with the 
following additions: 

1. NAMELIST/IN PT2/uses CECS72 variable names. 

2. hi the list below, the letters A, I, L, R refer to alphameric, integer, 


logical, real. 



Variable 

Type 

Description 

ADDI 

A 

' ADDF 

BAL 

L 

TRUE if equation is balanced 

BB 

A 

TRUE if equation is balanced 

CMOLE 

A 

* MOLE* 

CONV 

L 

TRUE if CECS72 matrix converges 

CL, CS 

L 

used in search for letter 

CS 

A 

* CAL 

DATA 

R 

temporary storage 

DATT 

R 

dummy variable 

DENOM 

R 

factor in extent of reaction calculation 

EPSL 

R 

extent of reaction 

FB 

A 

dummy variable 

FRAC 

R 

mole fraction 

GO 

L 

TRUE f CECS72 calculation is successful 

IC 

I 

dummy variable 

EDONE 

I 

message has been printed if = 0 

IEP 

I 

index in implied DO loop 

m 

I 

upper limit in implied DO loop 

n 

I 

reactant ratio 
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Variable 


Description 

m 

i 

number of temperature schedules 

IL 

i 

lower limit in implied DO loop 

ILAST 

i 

number of reactant ratios 

INS 

L 

TRUE if any products are solid or liquid 

IP 

I 

number of products 

m 

I 

index in implied DO loop 

IRR 

I 

number of .reactant ratios 

IS 

I 

index in DO loop 

IT 

I 

number of stops 

JL 

I 

number of temperature schedules 

JT 

I 

index in DO loop 

KB 

I 

base reactant 

KBB 

I 

not used 

KKLAST 

I 

number of products 

KKP 

I 

sum of number of products and reactants 

KP 

I 

number of products 

KR 

I 

index in DO loop 

LF, LL, LS 

I 

dummy variables 

LI, L2 

I 

dummy variables 

M 

I 

dummy variable 

N 

I 

dummy variable 

NCON 

J 

number of temperature schedules 

NOCON 

I 

dummy variable 

NOGO 

L 

TRUE if any step is unacceptable 

NPR 

I 

number of products 

NRE 

I 

number of reactants 

NRR 

I 

number of reactant ratios 

NT 

I 

reduced number of temperature schedules 
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Variable 

Type 

Description 

PMF 

R 

product mole fraction 

PR 

R 

pressure 

PRO 

A 

' PROD 

PROD 

R 

not used 

PROO 

R 

dummy variable 

REAC 

A 

'REAC* 

REACT 

R 

not used 

RMF 

R 

reactant mole fraction 

SB 

R 

temporary storage 

SIG 

R 

sum of mole fractions 

SING 

A 

* SIN* 

STOP 

A 

* STO» 

SI, S2, S3 

R 

temporary storage 

T, TB, TIB 

L 

used in character search 

TEM1, TEM2 

R 

temporary storage 

THM 

A 

i the» 

TT 

R 

dummy variable 

WARN 

A 

■ WAR' 

XDEN 

R 

factor in extent of reaction calculation 

XL 

L 

not used 

XLON 

R 

lowest acceptable extent of reaction 

XMUL, XNUM 

R 

factors in extent of reaction calculation 

XS 

A 

«s» 

XTRMAX 

R 

maximum allowed concentration 

XMF 

R 

mole fraction 
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1 



TABLE I. - DESCRIPTION OF CYCLES 



Reported (ref. 2) 


Revised 


L Mr(A) 







500 K 



600 K 


MfjO(s) + S0 2 + H 2 0 + S(f) 

m m 

Mj£0 4 <8) + H 2 S 

MgO(s) + S0 2 + HgO + S(f ) 


MfiS0 4 ( S > • n 2 s 


1400 K 



1400 K 


M*S0 4 (s) 



MgO(s) + SOg + 0. 5 0 2 

MgP0 4 (8) 


MfiO(s) i SO., t 0. 5 O., 


1300 K 



1300 £ 


lijS 


h 2 +S 

h 2 s 



ll 2 t 0. 5 S 2 





718 K 





0.5S 2 



m 

2. MtfB) 







F00 K 



600 K 


MgOfs) + S0 2 + H 2 0 + S(f ) 



M(^0 4 (s) + h 2 s 

MgO(s) + SO z +■ HgO + S(f ) 



MgS0 4 (s) + IlgS 


1400 K 



1400 K 


MbS0 4 (s) 

— 

MkO(s) + S0 2 + 0. 5 0 2 

MgS0 4 (s) 



MrO(s) + S0 9 + 0. 5 O., 


500 K 



500 K 


H 2 S + 0. 67 Bi(s) 

T ^ 

0. 33 Bi 2 S 3 (s) + H, 

HgS + 0. 67 BI(b) 



0. 33 Bi.,S.j(s) t II., 


1000 K 



1000 K 


0. 33 B£ 2 S 3 (s) 

m m 

0.67 Bi(f) + S 

0. 33 Bi 9 S 3 (s) 


0.67 10(1) i 0. 5 S„ 





718 K 





0.5S 2 

— 

S (1) 

3. Zn-Ba 







500 K 



500 K 


0. 25 BaS(s) + HgO 


0. 25 BaS0 4 (s) + H 2 

0. 25 BaS(s) + HgO 

— — 

0.25 B:iS0 4 (s) t 11., 


1400 K 



1400 K 


0.25 BaS0 4 (s) t0.5S 


0. 25 BaS(s) + 0. 5 SO z 

0. 25 BaS0 4 (s) + 0. 05 Sg + 0. 05 S g 


0.25 BaS(s) < 0. 5 SO., 


500 K 



500 K 


ZnO(8) +1.5 S0 2 


ZnS0 4 (s) + 0. 5 S(l) 

ZnO(s) i 1. 5 SO., 

— 

ZnSO.(s) +0. 5S(() 


1300 K 



1300 K 


ZnS0 4 (s) 

m 

ZnO(s) + SO z +0.50, 

ZnS0 4 (s) 


ZnO(s) + SO,, (0.5 <)., 





850 K 





0.5 S(f) 

— - 

0. 05 S., > 0. 05 S H 

4. Zn-Ca(A) 







500 K 



500 K 


0.25 CaS(s) + S(f) + H 2 0 


0. 25 CaS0 4 (s) + H 2 S 

0.25 CaS(s) + S(l) + HgO 

m m 

0. 25C:iSO,(s) ( II., S 


1400 K 



1400 K 


0. 25 CaS0 4 (s) + 0. 5 S 

: ^ ==: 

0. 25 CaS(s) + 0. 5 S0 2 . 

0. 25 CaS0 4 (s) + 0. 25 S g 


0. 25 CaS(s) ♦ 0. 5 SO., 


500 K 



500 K 


H.,S + 0. 67 Bl(s) 

= ^ ::=: 

0. 33 Bi 2 S 3 (s) + Hj, 

HgS + 0. 67 Bi(8) 


0. 33 Bi 2 S 3 (s) ( ll 2 


1000 K 



1000 K 


0. 33 K 2 S 3 (s) 


0.67 Bi(f) +S 

0. 33 BI 2 S 3 (8) 


0.67 Bi(t) ( 0. 5 S., 


500 K 



500 K 


ZnO(s) + 1. 5 S0 2 


ZnS0 4 (s) +0.5 S(l) 

ZnO(s) i 1,5 SOg 

— — 

ZnSO ( (s) • 0. 5 S (1) 


1300 K 



1300 K 


ZnS0 4 (s) 



ZnO<8) + SOg + 0. 5 Og 

ZnS0 4 (s) 


ZnO(s) i SOg • 0. 5 O., 





718 K 





0. 25 S 2 


0. 5 S(l) 
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tabu: I. 

- Continued. 



Reported (ref. 2) 


Revised 


5. Zn(A) 







850.K 



850 K 


IIjjO * 1. 5 S 


1I.,S t 0. 5 S0 2 

H 2 0 +0.55 S 2 + 0.05 S„ 

— ** 

l!.,S 0. 5 SO., 


1300 K 



1300 K 


It-jS 


H 2 1S • 

h 2 s 


II., * 0. 5S„ 


500 K 



500 K 


ZnO(s) + 1. 5 SO,, 

m m 

ZnSO j(s) + 0. 5S(I) 

ZnO(s) + 1.5 SO, 2 


ZnSOds) • 0. 5 S(f ) 


1300 K 



1300 K 


ZnS() 4 (s) 

M * 

ZnO(s) * SO., • 0. 5 O., 

ZnS0 4 (s) 



ZnO(s) ■ SO., • 0. 5 ()., 





850 K 





0. 5 S(i ) 



0.05 S., t 0.05 

6. Zn-Ca(B) 







500 K 



500 K 


0. 5 CiiS(s) + 2S({ > + 2II 2 0 


0. 5 CaSO,(s) + 2H 2 S 

0.5 CaS(s) + 2S(f ) t 2H.,0 

— - - 

0. 5 0aSO 4 (s) t 2I1.,S 


1400 K 



1400 K 


0. 5 CaS0 4 (s) + S 


0. 5 CaS(s) t- S0 2 

0. 5 CaS0 4 (s) + 0. 5 S., 

^ *“• 

0. 5 CaS(s) 4 SO., 


500 K 



500 K 


H 2 S 4 0. 67 Bi(s) 

— 

0. 33 Bi. 2 S 3 (s) + H 0 

H 2 S + 0. 67 Bi(s) 


0.33 Bi 2 S.,(s) i ll 2 


1000 K 



1000 K 


o. :t:i Bi 2 S. t (s) 

m w 

0.67 Biff) + S 

0. 33 Bi 2 S 3 (s) 


0.67 Bi(f ) + 0. 5 S., 


500 K 



500 K 


ZnS(s) + 2SO,, 

, * 

ZnS0 4 (s) i 2S(I) 

ZnS(s) + 2S0 2 


ZnS0 4 (s) t 2S(I) 


1300 K 



1300 K 


ZnS0 4 (s) 



ZnO(s) + S0 2 + 0. 5 0 2 

ZnS0 4 (s) 


Zn()(s) t SO,, • 0. 5 O., 


600 _K 



600 K 


ZnO(s) + »I,,S 



ZnS(s) + H 2 0 

ZnO{8) + H 2 S 



ZnS(s) i II., O 

7. Zn(B) 







850 K 



850 K 


2H 2 0 + 3S 

-i 9 

3 1 t.,S + St , 

2H z O + 0.75 S 2 + 0.2 S g 

— 

2II 2 S > S0 2 


1300 K 



1300 K 


h 2 s 


H 2 ‘S 

h 2 s 

-r m 

II., • 0. 5 S., 


500 K 



500 K 


ZnS(s) + 2S0 2 


ZnS0 4 (s) + 2S(() 

ZnS(s) + 2S0 2 

, 1 

ZnSO ,(s) f 2S(f ) 


1300 K 



1360 K 


ZuSOj(s) 


ZnO(s) * SOg + 0. 5 0 2 

ZnS0 4 (s) 

— 

ZnO(s) 4 SO,, i 0. 5 ()„ 


600 K 



600 K 


ZnO(s) + H 2 S 


ZnS(s) + H 2 0 

ZnO{s) + H 2 S 


ZnS(s) t ll.,<) 





850 K 





2S(J) 



0. 2 S 2 + 0. 2 S 8 

8. Fe-Ba 







500 K 



500 K 


0. 25 BaS(s) + H 2 0 

^ — . 

0. 25 BaS0 4 (s) + H, 

0. 25 BaS(s) + H 2 0 

m m 

0.25 BaS() 4 (s) < II , 


1400 K 



1400 K 


0. 25 BaS0 4 (8) + 0. 5 S 

m 

0. 25 BaS(s) t 0. 5 S0 2 

0. 25 BaS0 4 (s) + 0. 05 S 2 + 0. 05 S 8 

m 

0. 25 BaS(s) * 0. S( i ,, 


600 K 



400 K 


F»'.,t> 3 (8) * 2. 5 S0 2 


2FcS0 4 (s) t 0. 5 S(l ) 

Kc 2 0 3 (b) + 2. 5 S0 2 

^ 

2l eS(),(.s) * 0. 5 Sd ) 


1000 K 



1000 K 


2FeS0 4 (s) 


Fc,,0.,(s) + SO,, * SO., 

2FeS0 4 (8) 


l'e,,0.,|s) ' SO., SO., 


1300 K 



1300 K 


SO;, 

m m 

S0 2 + 0. 5 0 2 

so 3 

m 

SO,, < 0. 5 O, 





H30 K 





0. 5 S(l) 

- * 

0.05 S., ■ 0.05 K„ 


33 


TABLE L - Continued. 


Reported (ref. 2) 


9. Fc-(‘a 



0.25 CaS(s) + S<! ) + 11 2 0 

500 K 

0. 25 C 0 SO 4 (s) + 11,8 

0. 25 CaS0 4 <s> + 0. 5 S 

1400 K 

0. 25 CaS<s) t 0. 5 S0 2 

Fl',0.,(s) + 2. 5 SO., 

600 K 

2FeS0 4 (s) > 0. 5 S(l) 

2FcS0 4 (s) 

1000 K 

Fe 2 0 3 (s) + SO z + SOg 

H.,S >0.667 Bi(8) 

500_K 

0. 333 Bi^fs) + H 2 

0. 333 BigSgCs) 

1000 K 

0.667 Bl<f > +S 

SO., 

1300 

S0 2 + 0. 5 0 2 

10. Fe(A) 



H,0 + 1. 5 S 

850 K 

H 2 S *0.5 SOg 

H 2 S 

1300 K 

h 2 >s 

Kc.,0.,(s) + 2. 5 S0 2 

400 K 

2FeS0 4 (s) +0.5 S(l) 

2KoS0 4 (s) 

940 K 

Fe 2 0 3 (s) + SOg + SOg 

SO. ( 

1300 K 

SOg ♦ 0. 5 Og 

11. Fe(B) 



4H,,0 + 6S 

350 K 

4H 2 S + 2S0 2 

Fe.,0 3 (8) + 4H 2 S 

600 K 

2FeS 2 («) + 3H 2 0 + H 2 

2FcS 2 («) + 4S0 2 

500 K 

2FoS0 4 (s) + 6S(I ) 

2FoS0 4 

1000 K 

FCgOgfB) + SOg + SOg 

SO., 

1300 K 

SOg + 0.5 Og 


Revised 


0. 25 CaS(s) +S(f) + HgO 

500 K 

0. 25 t'a.S( i 

II., S 

0. 25 CaS0 4 (s) + 0. 25 Sg 

1400 K 

0. 25 CuS(s) * 0. 

5 SOg 

Fe 2 Og(B) + 2. 5 SOg 

400 K 

2I’oS0 4 (s) i 0. 5 

s<n 

2FeS0 4 (s) 

1000 K 

Kc 2 Og(s) i SOg 

* so., 

HgS +0.67 Bi(s) 

500 K 

0.33 Ui.,S.,(s) i 

Mg 

0.33 BigSgfS) 

1000 K 

0.67 Bi(l) * 0. 5 

S 2 

SOg 

1300 K 

SO., * 0. 5 O., 


0. 25 S 2 

718^ 

0.5 S(f) 



HgO + 0. 55 Sg + 0. 05 S 8 

850 K 

ll„.S t 0 5 SO, 

h 2 s 

l300 r K 

II, • 0. 5 S, 

Fe 2 Og(s) + 2. 5 SO, 

400 K 

2FeS0 4 (s) +0. 5S(») 

2FeS0 4 (s) 

1000 K 

Fe,0.,(8) * SO, « SO., 

SOg 

1300 K 

SO, + 0. 5 Og 

0. 5S(I) 

850 K 

0.05 S, * 0.05 S 8 

4HgO +0.6 S 2 +0.6 S 8 

850 

4llgS > 2SOg 

Fe 2 0 3 (8) + 4 HgS 

600 K 

2FcS,(s) + 3H.,() > II, 

2Fefi 2 (8) + 4SOg 

500 K 

2FeS() 4 (s) * 6S(0 

2FoS0 4 (8) 

1000 K 

Fc,0.,(s) i SO., • SO, 

SOg 

1300 K 

SOg >0.5 O, 

6S(I) 

850 K 

0.6 S, * 0.6 S 

•L n 


TABLE I. - Concluded. 
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TABLE IV. - COMPARISON OF EFFICIENCIES CALCULATED 
USING EQUILIBRIUM EXTENTS OF REACTION AND 


ASSUMING COMPLETE REACTION 


1 Process 

i 

i 

1 

Equilibrium 

Complete 
react! cm 

Equilibrium 



Complete 

reaction 

i 



1. Mg<A) 

0.50 

0.70 

0. 37 


3. Zn-Ba 

.71 

.80 

.46 

■ 1 

5. Zn(A) 

.52 

.66 


.40 

7. Zn(B) 

.45 

.62 

.26 

.38 

8. Fe-Ba 

.66 

.72 

.34 

.37 

10. Fe(A) 

.48 

.61 

.24 

.31 

! 11. Fe<B) 

.38 

.53 

.20 

.28 


TABLE V. - COMPARISON OF CALCULATED 


AND ESTIMATED EFFICIENCIES 
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